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1 . 0 INTRODUCTION 


This  report  summarizes  the  theory  and  operation  of  the  MCCABE  (Multi-Conductor-CABIE 
code)  which  is  an  analytical  tool  for  predicting  the  response  of  shielded  cables  or  cable 
bundles  in  X-ray  environments,  and  is  a multi-conductor  generalization  of  the  PICS  code 
reported  in  Ref.  9.  The  output  of  the  MCCABE  code  is  the  per-unit-leng'/n  eqivalent  circuit 
parameters  for  transmission  lines,  chief  among  which  is  the  radiation  driven  current  source 
terms.  The  underlying  assumptions  and  limitations  of  the  code  are  discussed  in  Sect. 2, 
but  the  main  assumptions  are  that: 

1)  the  cables  are  in  vacuum; 

2)  the  incident  flux  is  low  enough  such  that  radiation  induced  conductivity  and 
other  limiting  effects  may  be  ignored; 

3)  there  is  no  significant  x-ray  spectral  content  above  300  keV. 

The  MCCABE  formalism  presented  in  Sect.  2 differs  from  that  in  Ref.  1 in  that  the 

electron  transport  routines  hcve  been  completely  revised  to  use  the  analytic  charge 

4 

deposition  profile  analysis  recently  given  by  Dellin  and  MacCallum.  The  advantage  of 
their  deposition  profile  results  (besides  computing  speed)  is  that  electron  production 
from  both  materials  which  form  an  interface  (e.g.  conductor  and  dielectric)  is  automat- 
ically taken  into  account.  The  MCCABE  code  can  now  calculate  the  response  of  a cable 
made  from  materials  with  comparable  emission  efficiencies  (e.  . aluminum  and  Kel-F). 

The  operation  of  the  code  is  discussed  in  Sect.  3 where  sample  input  and  output  are 
also  presented.  One  important  feature  of  the  code  is  an  executive  subprogram  which  inter- 
rogates the  user  concerning  the  input  to  the  code  (cable  dimensions,  materials,  etc.). 

This  subprogram  sets  up  all  the  input  files  for  the  problem,  and  the  user  doesn't  have  to 
concern  himself  about  format.  The  running  time  for  a coax  problem  is  less  than  12  sec 
CPU  time  on  the  CDC  6600,  with  multi-conductor  cables  taking  longer. 

2) 

The  MCCABE  code  has  been  used  recently  In  a parameter  study  to  determine  the  sen- 
sitivity of  cable  response  to  various  cable  and  photon  source  parameters.  Before  running 
the  MCCABE  code  we  suggest  that  the  user  consult  Ref.  2 in  order  to  get  a rough  idea  of 
how  his  cable  will  respond. 

At  this  point  the  MCCABE  code  has  received  preliminary  experimental  verification  in 

7 8 

X-radiation  experiments  conducted  at  the  Simulation  Physics  facilities  ’ and  the  Aero- 

7 9 

space  Dense  Plasma  Focus  facility  ’ . References  7-9  concerned  either  coaxial  cable  re- 
sponse, or  common-mode  multi-conductor  cable  response;  Reference  1 considered  individual 

wire  response  of  multi-conductor  cable  response  as  well.  The  fluence  at  both  the  SPI  and 

-4  _3  2 

DPF  machines  was  on  the  order  of  10  to  10  cal/cm  . An  experimental  program  for  testing 
cables  is  presently  underway  at  the  SPI  facilities,  and  will  be  reported  later. 


As  mentioned  above,  the  code  is  applicable  to  the  region  where  radiation  induced 
conductivity  and  other  limiting  effects  are  unimportant.  We  believe  that  the  cross-over 
point  occurs  at  a few  tenths  of  a cal/cm  , and  that  beyond  this  the  response  is  sublinear. 
This  point  is  discussed  in  more  detail  in  Ref.  3.  The  implication  is  that  MCCABE  can 
still  provide  worst-case  estimates  of  cable  response  even  at  high  fluences. 


2.0  THEORY  OF  THE  MCCABE  CODE 


r 


r 


i 


l 

i 


The  problem  of  determining  the  response  of  shielded  cables  to  x-rays  can  be  consid- 
ered in  three  steps: 

Step  1:  Calculation  of  the  photon  transport,  electron  production  and  transport 

for  the  geometry  and  materials  of  the  cable  to  determine  the  profile  of 
deposited  charge  in  cable  dielectrics; 

Step  2:  Calculation  of  the  radiation  driven  source  terms  (Norton  equivalent 

drivers)  to  be  used  in  the  transmission  line  equations. 

Step  3:  Solution  of  the  transmission  line  equations  to  obtain  load  response. 

The  MCCABE  code  carries  out  the  first  two  steps,  again  assuming  that  no  limiting  occurs 
and  that  the  radiation  transport  is  independent  of  any  voltages  developed  along  the 
transmission  line.  With  these  assumptions  the  Norton  driver  can  be  inserted  into  trans- 
mission line  equations  in  a separate  step.  On  the  other  hand,  if  limiting  effects  are 
important,  then  steps  2 and  3 are  coupled  and  must  be  solved  simultaneously. 

In  the  subsections  below  we  discuss  the  following  topics: 

o transmission  line  equations  and  equivalent  circuit  modeling  (Sect.  2.1) 

o cable  geometry  and  photon  attenutation  (Sect.  2.2) 

o electron  deposition  in  cables  (Sect.  2.3) 

o the  solution  of  Laplace's  equation  in  two  dimensions  (Sect.  2.4) 
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2.1  TRANSMISSION  LINE  EQUATIONS  AND  EQUIVALENT  CIRCUIT  MODEL 

We  derive  the  canonical  transmission  line  equations  which  include  the  desired  X-ray 
induced  source  terms.  The  key  assumption  of  transmission  line  theory  is  that  only  the 
TEM  mode  of  propagation  is  important,  and  this  implies  that  the  electromagnetic  fields 
are  normal  to  the  axis  of  propagation.  An  immediate  consequence  of  this  is  that  the 
electric  field  is  derivable  from  a scalar  potential. 


Differential  Voltage  Transmission  Line  Equation: 

Consider  N wires  inside  a shield,  as  in  Fig.  2.1.1.  Consider  the  magnetic  flux 
linked  between  the  i^  conductor  and  the  chield  which  is  taken  as  reference.  From 
Lenz ' s law  we  have 


dA  . 


(2.1.1) 


The  contour  integration  is  defined  in  Fig.  2.1.1.  If  J,  the  photo  Compton  current 
density,  is  uniform,  and  normal  to  the  conductor  axes,  the  net  flux  linked  by  J is 
zero,  and  only  the  propagation  current  along  the  wires,  1^,  contributes  to  the  flux. 

By  definition,  each  of  these  will  make  a contribution  to  the  magnetic  flux  by  an 

amount  L..I.Az,  where  L,.  is  the  mutual  inductance  per  unit  length  of  the  wires  i and 

1J  J 1J  th  - 

j,  and  L is  the  self-inductance  per  unit  length  of  the  i wire.  Since  E is  normal 

to  the  conductor  axis  and  is  also  derivable  from  a scalar  potential,  the  LHS  is  just 

svi 

vi  + 4vi  - vi  : ST" 12  • <2-1-2> 


where  V^  is  the  potential  of  the  i wire  with  respect  to  reference,  and  we  obtain  the 
first  transmission  line  equation, 


3Vt  31. 

3z  ^ ^ij  3t 

j 


(2.1.3) 


Differential  Current  Transmission  Line  Equation: 

From  potential  theory  the  total  charge  per  unit  length  Q^  on  conductor  i is  given 
by 


total  charge 
per  unit  length 


d r q(r)i^i(r) 


induced  charge 
per  unit  length 


N 


capacitive  charge 


per  unit  length 


(2.1.4) 


(This  expression  is  derived  in  Appendix  A) . Here  q is  the  volume  charge  density  at 
any  instant  of  time,  ip^(r)  is  a solution  of  Laplace's  equation  with  unit  potential  on 
conductor  i,  and  zero  on  the  rest.  The  k„  are  the  elements  of  Maxwell's  capacitance 
matrix,  and  are  obtained  by  placing  a unit  potential  on  each  conductor  j grounding 
the  rest,  solving  Laplace's  equation,  and  obtaining  the  induced  charge.  Then  = k . 

Now,  the  rate  at  which  charge  is  deposited  in  the  element  Az  on  conductor  i,  is 
simply  the  difference  between  the  incoming  current  I^(z)  and  the  outgoing  current 
1^  (z+Az),  or  simply  -(3I/3z)Az.  This  is  equal  to  the  sum  of  two  terms:  (1)  (3Q^/3t)Az, 
which  includes  both  induced  and  capacitive  charge  (Eq.  2.1. A),  and  also  (2)  the  rate  of 
direct  charge  arrival,  Az  dl^  J*n^,  where  n^  is  the  surface  normal  directed  into 
the  conductor.  The  result  of  this  is  the  second  transmission  line  equation 


31.  N 3V . 

sr ' Ki  h,  it1  • 


(2.1.5) 


where  the  source  term  (i.e.,  the  Norton  equivalent  current  driver)  is  given  by 


Ki  =/d2r  lihi (?)  -/  dliH 


(2.1.6) 


and  3q/3t  is  related  to  the  photo-Compton  current  via  the  continuity  equation 


il  = _ v 

3t 


(2.1.7) 


By  construction  J only  includes  the  photo-Compton  current,  though  in  principle  a 
conductivity  term  could  be  added.  Since  we  have  assumed  that  both  J and  3q/3t  follow 
the  radiation  pulse,  then  so  do  the  K^. 

The  transmission  line  equations  (2.1.3)  and  (2.1.5)  may  be  solved  directly,  of 
course,  but  if  conventional  circuit  analysis  codes  are  to  be  used,  then  the  equivalent 
circuit  lumped  elements  of  a section  of  transmission  line  of  length  Az  are  shown  in 
Fig.  2.1.2.  If  N is  the  total  number  of  inner  wires,  the  model  consists  of  N current 
drivers  KiAz,N(N+l)/2  capacitances  C„Az,  and  twice  that  many  self-and  mutual  induc- 
tances l/2L^jAz.  The  elements  of  Maxwell's  capacitance  matrix  k_  are  related  to 
empirical  capacitances  C„  in  Fig.  2.1.2  by  the  relations 


Figure  2.1.2  Equivalent  circuit  model  for  a multi-conductor  transmission  line 
including  radiation  induced  source  terms. 
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2.2  MULTICONDUCTOR  CABLE  GEOMETRY  AND  PHOTON  ATTENUATION 

Next,  we  discuss  the  details  of  the  incident  photon  spectrum.  Throughout  this  work 
we  are  assuming  that  there  is  no  time  delay  in  electron-photon  transport  and  that  the 
current  drivers,  which  are  the  output  of  this  code,  are  proportional  to  the  instantaneous 
photon  flux.  With  this  in  mind  the  user  specifies  the  instantaneous  photon  flux  F and 
an  energy  spectrum.  For  convenience  this  could  be  the  peak  flux  associated  with  a given 
pulse.  Then  when  subsequent  transmission  line  codes  are  driven  by  the  Norton  equivalent 
drivers 

Norton  driver  = peak  Norton  driver  x pulse  envelope 

where  the  "pulse  envelope"  is  just  the  pulse  waveform  normalized  to  unity  at  the  peak 
of  the  pulse.  In  addition  all  results  scale  with  flux,  so  any  convenient  unit  of  flux 
is  appropriate. 

The  source  may  have  either  an  arbitrary  unnormalized  photon  energy  distribution 
U(E)  or  a blackbody  distribution  of  the  form 


U(E)  i 


exp(E/kT)  - 1 


(2.2.1) 


Suppose  a distribution  such  as  that  shown  in  Fig.  2.2.1  were  given. 


Figure  2.2.1  Unnormalized  photon  energy  distribution. 


One  proceeds  to  discretize  the  spectrum  by  picking  out  photon  energies  (the  dashed 
lines  in  the  figure)  which  are  not  necessarily  uniformly  spaced.  Then  the  total  un- 
normalized energy  flux  would  be  approximated  by 


r N 

/ U(E)  dE  *53  U(E  ) AE  , 
J j=i  3 3 


(2.2.2) 


E2  - Ei  » J = 1 » 


AV  2 (Ej+i " • j = 2>  N “ 1 * 


(2.2.3) 


en  - en-i  ’ j = N- 


Since  F is  the  total  energy  flux,  then  the  energy  flux  spectrum  F(E^)  Fs 


FU (E . ) FU(E.) 


F(Ej)  /U (E)  dE  ~ N 


Yj  U(E.)  AE. 


1 J 


(2.2.4) 


fhe  number  flux  spectrum  N(E^)  Fs 


F(E  ) 

»(E.)  - 


(2.2.5) 


and  the  number  flux  of  photons  with  energies  centered  about  E^  is 


n(E  ) = N(E.)  AE.  . 

i i 3 


(2.2.6) 


The  incident  photon  number  spectrum  n(E^)  now  must  be  attenuated  through  a set  of 
planar  shields  in  front  of  the  cable,  and  then  through  the  cable  itself.  It  is  assumed 


that  exponential  attenuation  is  sufficient  to  describe  the  situation;  in  other  words, 
multiple  photon  scattering  (from  Compton  scattering)  is  not  important  for  X-ray  spectra 
of  interest  (below,  say,  300  keV)  and  neither  detailed  photon  transport  has  to  be  taken 
into  account,  nor  do  "buildup  factors"  have  to  be  included.  In  brief, 


• oiit  * ™ 

n(E^)  = n(E^)  exp[-y(Ej)  x density  x g.p.l.]  (2.2.7) 

where  y is  the  attenuation  coefficient  and  g.p.l.  stands  for  geometrical  path  length. 

First  we  discuss  planar  attenuation  and  its  geometry.  We  had  in  mind  the  fact  that 
cables  will  rarely  be  directly  exposed  to  radiation  but  will  be  shielded  somewhat  by  a 
side  of  a box  or  the  skin  of  the  system,  for  example,  .‘dding  attenuators  is  also  a way 
of  hardening  the  system.  For  convenience  the  plane  of  the  incident  radiation  is  normal 
to  the  absorbing  plane,  but  the  photon  direction  itself  need  not  be  parallel  to  the 
absorbing  plane's  normal.  This  is  shown  in  Fig.  2.2.2. 
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The  geometrical  path  length  (g.p.l.)  is  clearly 

g.p.l.  = t sec  <j>,  (2.2.8) 

where  t is  the  attenuator  thickness  and  <J>  is  the  angle  that  the  photon  direction  makes 
with  the  attenuator  normal  (0  < <p  < 90°). 

Next  consider  the  attenuation  through  the  cable  itself  (Fig.  2.2.3).  Let  us  pick 
the  origin  of  coordinates  as  the  center  of  the  shield,  and  the  shield  axis  as  the  z axis. 
The  photons  are  considered  to  be  incident  in  discrete  planar  sheets,  parallel  to  the  x-z 
plane,  and  labeled  by  the  coordinate  y,  the  perpendicular  distance  of  the  plane  to  the 
z axis.  The  photon’s  direction  in  this  coordinate  system  is 


n , = i cos  d>  — k sin  d>  , 

phot 


(2.2.9) 


where  <p  has  the  same  meaning  as  before  in  our  discussion  of  planar  attenuation.  What  is 
the  geometrical  path  length  of  a photon  starting  at  the  point  P^(x^,y^,z^)  and  arriving 
at  ?2^x2  = ^1,Z2^'  because  of  the  particular  choice  of  axes,  y0  = y , the  projection  of 
the  geometrical  path  length  onto  the  x axis  is  given  by 


g.p.l.  = projected  g.p.l.  x sec  . 


(2.2.10) 


This  will  be  true  for  any  sheet  of  photons  regardless  of  its  coordinate  y.  The  point  of 
this  is  that  we  can  1)  talk  about  projected  g.p.l. 's  exclusively,  as  if  the  radiation 
were  normally  incident  (along  the  x-axis),  and  2)  reduce  the  problem  to  two  dimensions. 
After  calculating  the  projected  g.p.l.,  we  simply  have  to  remember  to  multiply  by  sec 

So  now  we  look  at  a cross  section  of  a multiconductor  cable,  as  in  Fig.  2.2.4.  A 
plane  containing  the  incident  photons  intersects  the  cable  at  coordinate  y.  The  problem 
of  finding  the  projected  path  length  reduces  to  1)  finding  the  x-intersections  of  the 
line  y with  the  different  interfaces  (the  dots  in  the  Fig.  2.2.4)  and  2)  taking  the  dif- 
ference between  successive  intersections.  In  terms  of  the  notation  in  the  figure,  these 
x-intersections  for  a fixed  y are  given  by 


VPF 


x = p cos  6 ± \l  l 0 \ - [y-P  sin  6] 


(2.2.11) 
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The  radius  a is  to  be  used  for  conductor  intersections,  b for  dielectric  intersections. 
(For  purposes  of  calculating  attenuation,  we  ignore  gaps).  If  an  intersection  is  to 
occur,  it  occurs  in  pairs.  If  the  discriminant  is  negative,  no  intersection  occurs.  In 
summary,  we  calculate  x-intersections  via  Eq.  (2.2.11),  calculate  the  projected  path 
length  by  taking  differences,  multiply  by  sec  <j>  to  set  the  full  geometrical  path  length, 
"remember"  what  material  it  is  in,  look  up  the  appropriate  attenuation  coefficient,  and 
attenuate  the  number  flux  n(Ej)  to  that  point  via  Eq.  (2.2.7). 

For  purposes  of  calculating  electron  deposition  in  the  dielectric  it  is  necessary  to 
know  the  direction  cosines  of  the  normal  to  the  conductor  at  the  point  where  the  plane 
(labeled  by  y)  intersects  the  conductor  at  x.  Assuming  the  intersection  r 5 (x,y)  is 
known,  the  conductor  normal  unit  vector  satisfies  (see  Fig.  2.2.4), 


p _ i(x  - p cos  9)  + j (y  - p sin  6) 


iy  + jy 

x J y 


(2.2.12) 


which  identifies  the  direction  cosines  y and  y for  interior  conductors.  For  the  shield 

x y 

the  direction  cosines  are  y = -x/a,  y = -y/a  respectively.  The  angle  that  the  incident 

x y 

photon  makes  with  the  conductor  normal  is 


cos 


phot 


n = cos 


-1 


(y  cos  4> ) 


(2.2.13) 


where  Eqs.  (2.2.9)  and  (2.2.12)  have  been  used. 


2.3  ELECTRON  DEPOSITION 

To  evaluate  Eq.  (2.1.6)  for  the  Norton  drivers  we  require  the  emission  current  from 

each  conductor  surface  and  the  deposited  charge  profile  in  the  dielectric  volume.  To 

4) 

obtain  this  we  apply  the  analytic  prescriptions  given  by  Dellin  and  MacCallum  for  the 
case  of  a planar  interface  between  two  semi-infinite  media.  While  the  conductor- 
dielectric  interfaces  in  real  cables  are  not  planar,  the  distance  traveled  by  electrons 
is  usually  much  smaller  than  conductor  dimensions,  and  the  planar  approximation  is 
reasonable. 

The  first  term  in  the  expression  for  the  Norton  driver  (Eq.  2.1.6)  is  the  charge 
deposition  term  involving  3q/3t.  Since  charge  deposition  occurs  in  the  region  surround- 
ing each  conductor  and  does  not  extend  very  far  into  the  bulk  dielectric,  at  least  for 
x-rays  whose  energy  is  less  than  3 00  keV,  it  is  convenient  to  break  up  the  volume  integra- 
tion over  the  entire  dielectric  into  a set  of  integrations  near  each  conductor  surface. 


/d2r  S *i  = £ / d\  /dn  ^ *i 


(2.3.1) 


Here  di^  is  a element  of  arc  length  about  conductor  k and  n is  integrated  along  the 
conductor  normal  from  the  conductor  - dielectric  interface  into  the  volume.  Dellin  and 
MacCallum  have  found  that  the  deposition  profile  exponentially  decays  away  from  the 
interface  value,  i.e. 


H ^ Q,‘0)2 

. y >1 

dt  h ;<;> 


' (0) 
Qfi  n 


(2.3.2) 


•(0)  * (1) 

All  quantities  in  Q;  and  Q),  are  obtained  from  the  QUICKE2  code.  However,  the  emission 
N F B ^ 

efficienceis  J ,J  , and  J which  in  QUICKE2  have  the  units  electrons/photon,  are  multiplied 

ii..  * (0) 

by  -|e|n,  where  n is  the  attenuated  number  flux  of  photons  in  Eq.  (2.2.7),  and  Q has 

2 1J  2 

the  units  of  amps/cm  . Also,  the  electron  ranges  in  QUICKE2  which  have  the  units  g/cm 

are  converted  to  cm,  and  Q^d^has  the  units  of  amps/cm. 

The  second  term  in  Eq.  (2.1.6)  is  evaluated  as 


/ h • -fii  i j; 


(2.3.3) 
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where  J is  the  normal  component  of  the  interface  current.  (The  change  of  sign  comes 

1 X 

from  the  fact  that  is  directed  out  of  the  conductor  into  the  dielectric  in  contrast 
to  n^).  is  easily  obtained:  the  total  charge  deposition  rate  is  » but 

this  must  be  equal  to  the  difference  between  current  coming  into  the  volume,  i.e.,  the 

I N 

interface  current  J^,  and  the  current  leaving  it,  i.e.  the  net  bulk  current,  J^.  We 

obtain 


j]  =£  Q<0)  + J? 


1 jTi 


(2.3.4) 


,N  . 


where  is  taken  from  QUICKE2  data. 


The  above  expressions  apply  for  the  case  of  mono-energetic  incident  X-rays.  For 
the  case  of  a distributed  source,  the  full  Norton  drivers  must  be  built  up  by  super- 
position. The  first  index  i = 1,  2 indicates  the  medium  in  which  the  deposition  occurs 

and  the  index  1,  2 indicates  the  medium  in  which  the  electrons  are  created.  The  numbers 

* (0) 

1 and  2 refer  to  the  media  traversed  first  and  second  by  the  X-rays.  Q..  represents  the 

"(1) 

total  charge  deposition  rate  in  the  given  region,  while  Q . is  the  first  moment  of  it* 


Dellin  and  MacCallum  give  the  following  expressions  for 


Q<?>  and  Q<}>: 
ij  il 


*(0)=  -(-1)^  cosdyJT.  - (-1)1+^  (J^  + J®)  yjr  r2yj3/2  , (2.3. 

i 1 ^ * 1 

+ yp~2 


3) 


*(D_  1 riRi(Ej) 


V 1+1/2  S3_. 


\l  1+1/2  S,  + V 1+1/2  S„ 
* 1 ^ 


(-1)1R  (E . ) JN  coseVI/2 

. J - ,J  .J ■ - 


V 1+1 / 2 Sj 


-51J  + 


(2.3.4) 


*In  general,  there  are  sets  of  Q^^and  Q^^for  each  electron  source  type  (photo-electric. 
Auger,  Compton),  but  because  of  our  restruction  to  photon  energies  below  300keV  where 
the  photo-electric  effect  dominates,  we  calculate  only  one  set  of  the  and 
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2.4  LAPLACE  EQUATION  IN  TWO  DIMENSIONS 

We  present  a technique  for  solving  the  two-dimensional  Laplace  equation  which  takes 
advantage  of  the  fact  that  a cross-section  of  a multiconductor  cable  shows  a collection 
of  arbitrarily  placed  circles,  as  in  Fig.  2.2.4.  Rather  than  go  through  the  laborious 
procedure  of  finite  differencing  the  Laplace  equation,  we  simply  write  down  the  formal 
solution, 


11  „ 

4tte  (r ) = >7  a.  /*  11  dQ'  0.(9')  £n 

i=l  XJo 


(2.4.1) 


where 


r'  = P±  + a',  and  p1  = (pi>  6J,  a'  = (a.,,  0'), 


The  various  geometrical  quantities  are  shown  in  Fig.  2.2.4.  But,  in  this  development, 
we  assume  that  the  dielectric  is  homogeneous  and  that  there  are  no  gaps.  The  logarith- 
mic term  is  seen  to  be  the  Green's  function  for  a cylinder  of  radius  a^  which  is 
infinite  in  length:  this  guarantees  not  only  that  the  potential  satisfies  Laplace's 
equation,  but  that  it  also  vanishes  on  the  shield  whose  radius  is  a^.  The  unknown 
charge  density  on  each  conductor  includes  both  free  and  bound  charge  (free  charge 
density  = Ko^  where  K is  the  dielectric  constant  near  conductor  i) . Since  we  are  assum- 
ing a homogeneous  dielectric,  there  are  no  dielectric-dielectric  interface  charge  den- 
sities to  worry  about.  (The  complication  of  dielectric  interfaces  is  discussed  below). 
The  first  step  in  obtaining  the  unknown  charge  densities  is  to  expand  Eq.  (2.4.1)  in 
circular  harmonics  (cos  n0',  sin  n9 ' ) ; 


N 00 

4ire  <K?)  £ V (r)-X  , (2.4.2) 

o . , in  in 

i=l  n=o 


where  the  spinor  is  given  in  Appendix  B.  The  spinor  has  replaced  the  unknown 
charge  densities  by  a set  of  unknown  coefficients, 


X 


in 


2tt 

de'o  (e') 


(2.4.3) 


Che  very  first  one  of  which  yields  exactly  the  induced  free  charge  per  unit  length  = 


K X 


io 


The  coefficients  X are  determined  by  placing  the  prescribed  potentials 


(p  + a^) , a^  = (a.,  8),  on  the  boundaries,  and  multiplying  Eq.  (2.4.2)  by 


(4ir) 


-1/. 


2n 


d8  (cos  n0,  sin  n0)  . This  results  in  the  matrix  equation 


£ 


Vn* 


(2.4.4) 


where  the  matrix  A is  given  in  Appendix  B.  (In  practice  the  infinite  sum  over  n'  is 
truncated).  This  procedure  determines  the  capacitance  matrix  as  well. 


Correction  for  Non-Homogeneous  Dielectric: 

Multiconductor  cables  in  practice  will  not  usually  have  homogeneous  dielectrics; 
contrary  to  our  assumption  above,  the  existence  of  gaps  between  insulation  and  conductors 
as  well  as  dielectric-dielectric  interfaces  is  the  more  likely  possibility.  Our  assump- 
tion is  that  non-homogeneities  in  the  dielectric  will  not  affect  the  capacitance  matrix 
very  much  (<50%).  On  the  other  hand  the  effect  on  the  potential  distribution  near  con- 
ductor-dielectric interfaces  is  considerable,  as  will  be  made  clear  by  the  following 
argument.  Consider  a parallel  plate  capacitor  whose  plate  separation  is  D,  but  with  a 
small  gap  d (d<<D)  which  nonetheless  is  larger  than  the  range  p of  an  electron  penetrat- 
ing the  dielectric.  If  the  dielectric  medium  were  uniform  then  the  potential  distribu- 
tion at  X = d + p~d,  would  be  4> /V  = d/D.  However,  the  presence  of  a gap  (whose  dielec- 
tric constant  was  unity)  would  imply  a potential  at  the  same  point  of  4>/V  = Kd/D,  again 
assuming  d<<D.  This  suggests  that  for  electrons  which  cross  gaps,  we  calculate  the 
Laplace  solution  using  the  formalism  of  this  section  where  the  dielectric  constant  is 
that  of  the  dielectric  insulation,  and  then  scale  the  charge  profile  by  multiplying  by 
K.  (In  principle,  of  course  one  can  extend  the  formalism  of  the  last  section  by  demand- 
ing the  normal  component  of  the  displacement  vector  to  be  continuous  at  the  dielectric- 
dielectric  interface.  This  would  introduce  a set  of  dielectric  interface  charge  densities). 
The  i^(r)in  Eq.  (2.1.6)  are  obtained  by  setting  = 1 in  Eq.  (2.4.4)  and  then  using  Eq. 
(2.4.2)  to  obtain  the  potential. 


3.0  STRUCTURE  AND  USE  OF  THE  MCCABE  CODE 


The  MCCABE  code  was  written  for  use  on  an  interactive,  time-share  computer  system. 
The  code  itself  consists  of  four  separate  subprograms,  each  of  which  can  be  run  separate- 
ly. Briefly,  these  subprograms  accomplish  the  following: 

1.  FLUXPAK  concerns  itself  with  the  incident  x-ray  spectrum  and  looking  up  the  ap- 
propriate photon  attenuation  coefficients  as  a function  of  this  energy  spectrum. 

2.  ATTNPAK  calculates  the  photon  attenuation  of  the  incident  spectrum  until  it 
reaches  conductor/di electric  interfaces. 

3.  CIRCLES  calculates  the  capacitance  matrix  and  solves  Laplace's  equation. 

4.  DRIVPAK  calculates  the  electron  deposition  in  the  dielectric,  and  the  appropri- 
ate Laplace  equation  solutions  to  obtain  the  Norton  equivalent  drivers. 

Needless  to  say,  all  these  subprograms  are  communicating  with  each  other  in  a com- 
plicated way.  Fig.  3.1  gives  a flowchart  of  input  and  output  between  the  user  and  the 
various  subprograms.  Table  3.1  lists  the  information  contained  in  each  input/output 
file. 


Because  of  the  complicated  communication  between  subprograms  we  have  written  the 
MCHELP  code  which  handles  input  and  program  control  without  the  user  having  to  think  too 
much  about  it.  MCHELP  is  a user-oriented  code  which  interrogates  the  user  and  then  per- 
forms two  functions: 

1.  It  constructs  all  necessary  input  files  from  user-supplied  data  on  the  cable 
geometry  and  x-ray  environment. 

2.  It  sets  up  a job  control  file  for  running  the  various  subprograms. 

Fig.  3.2  gives  an  example  of  a session  between  user  and  MCHELP  for  the  case  of  a simple 
coax.  The  explanatory  text  provided  by  MCHELP  identifies  the  variables,  and  their  units, 
as  required  by  MCCABE.  The  job  file  (TAPE44)  assumes  that  the  various  input  files  (TAPE5, 
TAPE3,  TAPE15)  have  been  given  arbitrary  file  names  and  stored  in  memory.  Furthermore, 
it  is  assumed  that  relocatable  binary  versions  of  FLUXPAK,  ATTNPAK,  CIRCLES,  and  DRIVPAK 
have  been  stored  under  the  names  BFLUX,  BATTN,  BCIRCL,  and  BDRIV,  respectively.  Fig. 3. 3 
gives  the  job  file  produced  during  the  session  shown  in  Fig.  3.2.  The  output,  of  this 
MCCABE  run  is  shown  in  Fig.  3.4.  This  problem  required  12  CPU  seconds  on  a CDC  6600. 

Fig.  3.5  lists  the  input  files  for  a sample  three-wire  multiconductor  cable  as  produced 
by  MCHELP.  Fig.  3.6  gives  the  output  of  the  CIRCLES  and  DRIVPAK.  subprograms  of  MCCABE 
which  contains  the  calculated  capacitance  matrix  and  Norton  drivers.  This  problem  re- 
quired 87  CPU  seconds. 


Table  3.1  Data  files  used  with  MCCABE 


Name 

Use 

TAPE5 

TAPE3 

TAPE15 

TAPE6 

TAPE4 

TAP El 6 

TAPE24 

TAPE8=NSTAB 

TAPE9 

TAPE27=DELMAC 

TAPE10 

TAPE12 


input  of  X-ray  spectrum  and  of  cable  materials  description 
input  of  cable  geometry 

input  to  control  Laplace  equation  solution 

output  of  photon  cross  sections 

output  of  attenuated  photon  flux  at  interfaces 

output  of  checks  on  coefficients  of  the  potential  expansion 
and  capacitance  matrix  calculation 

output  of  electron  transport  and  Norton  driver  calculation 
photon  cross  section  data  for  FLUXPAK  (binary) 
photon  cross  sections  selected  for  ATTNPAK 
photo-Compton  current  data  from  QUICKE2  for  DR1VPAK 
photon  flux  at  conductor/dielectric  interface  for  DRIVPAK 
coefficients  of  potential  expansions  for  DRIVPAK 
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FLUXPAK 


TAPE6  TO  OUTPUT 


Figure  3.1 


MCHELP 

(SETS  UP 

INPUT 

FILES, 

TAPE5, 

TAPE3, 

TAPE  15) 


TAPE8=NSTAB 


W 


TAPE3 


I M 


TAPE15 


1 

f 

TAPE9 

1 

r 

ATTNPAK 

CIRCLES 


TAPE  12 


TAPE27=DELMAC 


TAPE4  TO  OUTPUT 


TAPE16  TO  OUTPUT 


TAPE10 


DRIVPAK 


TAPE26  TO  OUTPUT 


Flow-chart  of  input/output  communication  among  subroutines  in  MCCABE. 


MCHELP  8/3/7? 

THIS  CODE  SETS  UP  INPUT  INFO  FOR  THE  MCCABE  CODE 
THE  INFO  REQUIRED  IS  DIVIDED  INTO  2 PARTS: 

CABLE  GEOHETRY  AND  SPECTRUM 
EXECUTION  INSTRUCTIONS 

THE  OUTPUT  OF  THIS  CODE  IS  A PERFORH  FILE(TAPEAA) 
AND  VARIOUS  FILES  WHICH  NAY  BE  SAVED 


CABLE  GEOMETRY  AND  SPECTRUM 

INPUT  INFO  UHICH  YOU  GIVE  SHOULD  FOLLOU 
NAMELIST  RULES 

WITH  END  OF  ENTRY  INDICATED  BY  * 

IF  YOU  HAVE  INFO  IN  PERM  STORAGE 
PROVIDE  DUMMY  INFO  OR  NONE  AT  ALL 

DO  YOU  UANT  TO  CONSTRUCT  TAPE5 
FOR  FLUXPAK?  < 1 =YES,0=NO > 

? 1 


PHOTON!  NAMELIST 

FLUX=FLUX  IN  CAL/CH2SEC 

NNIX=  NUMBER  OF  DISTINCT  MATERIALS (MAX3! 0) 

N=  NUMBER  OF  PHOTON  ENERGIES  (NAX=SO) 

EPHOT (I ) ,1=1 ,N  ARE  PHOTON  ENERGIES(KEV) 

U(  I ) , 1 = 1 ,N  ARE  UNNORM.  ENERGY  SPECTRA(/KEV) 

IF  EQUALLY  SPACED  ENTER  EPHOT(l)  AND  EPHOT (N)  ONLY 
XKT  IS  BB  TEMP  IN  KEV(IGNORED  IF  NOT  ENTERED) 

PHI  IS  ANGLE(DEGREES)  BETWEEN  SHIELD  NORMAL 
PROJECTED  IN  PLANE  OF  INCIDENT  RADIATION 
AND  DIRECTION  OF  INCIDENT  RADIATION  <<90  DEG) 
IPRT3NONZERO  GIVES  LONG  PRINT 


BPH0T0N1 

? FLUX=1 ,N=1 ,EPHQT350fU=1 fNMIX=A* 

YOU  WILL  BE  ASKED  FOR  4 

MIXTURE  NAMELIST 

FOR  EACH  MATERIAL (OF  UP  TO  10  COMPONENTS): 

THE  ORDER  OF  ENTRY  DOES  NOT  NATTER 

BUT  THE  MATERIAL  LABEL  INDEX 

WHICH  FOLLOWS  LATER  IN  THE  INPUT  CORRESPONDS 

TO  THIS  ORDER  OF  ENTRY 

ZA=  Z NUMBERS! IN  ORDER  OF  INCREASING  Z) 

WT3  WEIGHT  FRACTIONS 
ATMTONIC  WEIGHT 
DENSITY=DENSITY<G/CM3 ) 

IZDM=Z  NO.  OF  ELEMENT  OR  LABEL  NO.  IN  DELMAC  TABLE 


Note: 


Dielectric  type 

Label  No. 

water 

901 

teflon 

902 

PVC 

903 

Mylar 

904 

Kapton 

903 

Halar 

906 

Kel-F 

907 

Fieure  3.2  Terminal  printout  from  MCHELP  for  a simple  coax. 
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♦MIXTURE 

? ZA=29,AT=43.54,UT=1  ,DENSI  TY=8.9 , IZDM=2<5» 

♦MIXTURE 

? ZA=4,9, AT=12 .01 , 19. ,UT= .24, . 74, DENSITY=2.2, IZDM=9Q2$ 
♦MIXTURE 

? ZA=50,AT=11S.7,UT=1 ,DENSITY=7.2 , IZDN=50* 

♦MIXTURE 

? ZA=47,AT=107.9,UT=1,DENSITY=10.5,IZDM=47« 


TAPE5  HAS  BEEN  CONSTRUCTED 


DO  YOU  UANT  TO  CONSTRUCT  TAPE3  AND  TAPE15 
FOR  ATTNPAK  AND  CIRCLES’  ( 1 =YES, 0=N0 ) 
f 1 


ATTPLN  NAMELIST 

NPLANE-  NO.  OF  ATTENUATION  PLANE S ( MA Y BE  0) 

T( I ) , 1=1 , NPLANE  = PLANE  THICKNESS  IN  CM 
IPLANLB,  1 = 1 ,NPLAN£=MIXTURE  INDICIES 
UHICH  IS  THE  ORDER  IN  UHICH  MIXTURES  WERE  READ 
♦ATTPLN 
? i 


CONDUCT  NAMELIST 

ALL  DISTANCES  IN  CM,  ANGLES  IN  DEGREES 
NCON=NO.  OF  COND'S  INCLUDING  SHIELD 
TSHIELD=SHIELD  THICKNESS 
RADC ( I ) , 1 = 1 , NCON  ARE  COND.  RADII 
RADD  ARE  DIELECTRIC  RADII 

GAP  ARE  GAP  SIZES 

THE  FIRST  ENTRY  OF  RADC, RADD, GAP  REFERS  TO  SHIELD 
MATLBLC(I), 1=1, NCON  IS  THE  MATERIAL  LABEL  INDEX 
MATLBLD  < I ) 

MAT L BL F ( I ) 

UHERE  C IS  THE  CONDUCTOR,  D DIEL . , F FLASHING, 

AND  THE  LABEL  INDEX  CORRESPONDS  TO 

THE  ORDER  IN  UHICH  MIXTURES  UERE  READ  IN  ABOVE 
MATBAK  IS  THE  BACKGROUND  DIELECTRIC  INDEX 
DFLASH(I), 1 = 1 , NCON  ARE  FLASHING  THICKNESSES 
RO(I),THET(I), 1=1, NCON  ARE  THE  COORDS.  OF  COND.  I 
2 TIMES  NY  IS  THE  NO.  OF  Y PLANES  INTERSECTING  SHIELD 
IPRT=NONZERO  GIVES  LONG  PRINT 


♦CONDUCT 

? NC0N=2,R0=0,0,NY=10 
? MATLBLC=1 ,1 ,MATL8LD=2,2,NATBAK=2 
? RADC=. 1,-0299 
? GAP=.005, .001 
’ RADD=. 0451,. 045 
? TSHIELD=.01 

’ MATLBLF=3,4,DFLASH=.0002, .00021 


TAPE3  HAS  BEEN  CONSTRUCTED 


INPUT  TO  CIRCLES  AND  TAPE  1 5 CONSTRUCTION 

MOST  INPUT  DAS  GIVEN  IN  ATTNPAK  FOR  CIRCLES 

PERM  NAMELIST 

PER=DIELECTRIC  CONSTANT 

MMAX=MAX  TERMS  IN  HARMONIC  EXPANSION 

EPSLN  TIMES  MIN.  DIAG  ELEMENT  IS  LARGEST  OFF-DIAG 

♦ PERM 

? PER=2.1,MMAX=1* 

TAPE15  HAS  BEEN  CONSTRUCTED 


/ 


Note : 


Note : 


MMAX=1  is  sufficient 
for  coaxes.  MMAX=4-6 
is  sufficient  for  cable 
bundles  depending  on 
degree  of  assymetry. 
EPSLN=0  means  no  off- 
diagonal  matrix  ele- 
ments are  discarded 
before  entering  sparse 
matrix  routines 


The  flashing  modification 
is  a phenomenological  ad- 
justment of  the  interface 
current  discussed  in  Refs. 
2 and  7. 


Figure  3.2  (continued)  Terminal  printout  from  MCHELP  for  a simple  coax. 
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EXECUTION  INSTRUCTIONS 


SHOULD  FLUXPAK  BE  EXECUTED 
( 1=YES,0=N0) 

? t 

NAME  OF  TAPE5 

? C0AX5  

NAHE  OF  TAPE9 
? C0AX9 


_ Note:  These  file  names 
are  arbitrary. 


SHOULD  ATTNPAK  BE  EXECUTED 
< t=YES,0=N0) 

? 1 

NAHE  OF  TAPE 3 

? C0AX3  

NAHE  OF  TAPE9 

? C0AX9  

NAHE  OF  TAPE10 
? COAX  10  


SHOULD  CIRCLES  BE  EXECUTED 
(1=YES,0=N0) 

? 1 

NAHE  OF  TAPE  15 

? COAX  15  

NAHE  OF  TAPE12 
? C0AX12  


SHOULD  DRIVPAK  BE  EXECUTED 
( 1 = YES  r 0=N0 ) 

? 1 

NAHE  OF  TAPE10 

? COAX 10  

NAHE  OF  TAPE  12 
? COAX  12  


FILES  CONSTRUCTED: 

TAPE441PERF0RM  FILE) 

TAPE51  INPUT  FOR  FLUXPAK) 

TAPE3( INPUT  FOR  ATTNPAK) 

T APE  1 5 < INPUT  FOR  CIRCLES) 

LIST  TAPE44  TO  SEE  NAMES  ASSIGNED  TO  DATA  FILES 
REPLACE  TAPE5, TAPE3, TAPE  1 5 UNDER  THIER  NEU  NAHES 
THEN:  PERFORM, TAPE44 


Figure  3.2  (continued)  Terminal  printout  from  MCHELP  for  a simple  coax. 
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GET ,TAPE5=C0AX5 
GET ,TAPE8=NSTAB 
GET,BFLUX 
BFLUX,TAPE4=0UTPUT 
REPLACE, TAPE9=C0AX9 
GET, TAPE3=C0AX3 
GET , TAPE9=CQAX9 
GET,BATTN 
BATTN,TAPE4=0UTPUT 
REPLACE, TAPE 10=C0AX1 0 
GET,TAPE15=C0AX15 
GET  ,BCIRCL 

BCIRCL,TAPE16=0UTPUT 
REPLACE, TAPE12=C0AX12 
GET,TAPE10sC0AXf 0 
GET ,T  APE  1 2=C0AX 1 2 
GET,TAPE27=DELNAC 
GET  ,BDRIV 

BDRIV,TAPE26=0UTPUT 


Figure  3.3 


Job  file  created  by  MCCABE  for  the  example  given  in  Figure  3.2. 
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* _tl 


1  FLUXPAK  8/3/77 


NO.  OF  MIXTURES 
FLUX ( CAL/CM2SEC ) 

PHI ( DEGREES) 

PHOTONS/CAL 
TOTAL  NFLUX < /CM2SEC ) 
INDEX/NUMBER  FLUX(/CM2SEC) 
1 5.23E  + 1 4 


4 

I.OOOOOE+OO 

0. 

5 . 22600E+ 1 4 
5.22600E+1 4 


MIXTURE  NO./Z  NO. /DENSITY/ATONIC  UT. 

1 2.90E+01  8.90E+00  6.35E+01 

MIXTURE  NO./Z  NO. /DENSITY/ATOMIC  UT. 

2 8.28E+00  2.20E+00  1.73E+01 

MIXTURE  NO./Z  NO. /DENSITY/ATOMIC  UT. 

3 5.00E+01  7.20E+00  1.19E+02 

MIXTURE  NO./Z  NO. /DENSITY/ATOMIC  UT. 

4 4.70E+01  1.05E+01  I.08E+02 


1 DRIVPAK  8/3/77 


GEOMETRY  VARIABLES  TRANSFERRED  FROM  ATTNPAK 

NCON*  2 

RADC*  1.00E-01  2.99E-02 

GAP=  5.00E-03  1.00E-03 

RADD*  6.51 E -02  6.50E-02 

R0=  0.  1.00E-07 

THET-  0.  1 0. 

MATLBLC=  1 I 

MATLBLB=  2 2 

MATLBLF*  3 4 

DFLASH*  2.00E-04  2.00E-04 

PHOTON  ENERGIES 
NE-  1 

EPHOT=  5.00E+01 

TOTAL  EMITTED  CURRENT  FROM  SHIELD  AND  EACH  CONDUCTOR 
I SURF  = 8.36E-08  2.40E-08 

TOTAL  INDUCED  CURRENT  ON  EACH  CONDUCTOR  DUE  TO  DEPOSITED  CHARGE 
I AREA=  -2.97E-08 

NORTON  DRIVER  TO  EACH  CONDUCTOR! SHORT  CIRCUIT  CURRENT)  (A/CM) 
ISC*  -5.67E-09 


Figure  3.4 


MCCABE  output  for  the  coax  defined  in  Figure  3.2 


1 ATTNPAK  8/3/77 


X/Y/X  DIR/YDIR/ANGLE/CON  NO. 

ENERGY  RIN/NO.  FLUX  UNIT  LENEGTH 

-3.12E-02  -9.50E-02  3.12E-01  9.50E-01  1.25E+00 

1 2.7BE+1 3 

3.12E-02  -9.50E-02  -3.12E-01  9.50E-01  1 .89E+00 

1 2.71E+13 

-5.27E-02  -8.50E-02  5.27E-01  8.50E-01  1.02E+00 

1 1.38E+13 

5.27E-02  -8.50E-02  -5.27E-01  S.50E-01  2.13E+00 

I 1.33E+13 

-6.61E-02  -7.50E-02  6. 61 E-01  7.50E-01  8.48E-01 

1 1.1AE+13 

A.A1E-02  -7.50E-02  -6.6tE-0t  7.50E-01  2.29E+00 

1 1 . tOE+13 

-7.60E-02  -A.50E-02  7.60E-01  A.50E-01  7.08E-01 

I I.04E+I3 

7.A0E-02  -6.50E-02  -7.60E-01  6.50E-01  2.43E+00 

1 9.81E+12 

-8.35E-02  -5.50E-02  8.35E-01  5.50E-01  5. 82E-0 1 

1 9.70E+12 

8.35E-02  -5.50E-02  -8.35E-01  5.50E-01  2.56E+00 

I 9.0AE+12 

-8.93E-02  -4.50E-02  8.93E-01  4.50E-01  4.67E-01 

1 9.20E+12 

8.93E-02  -4.50E-02  -8.93E-01  4.50E-01  2.A7E+00 

1 8.55E+1 2 

-9.37E-02  -3. 50E-02  9.37E-01  3.50E-01  3.58E-01 

1 8.86E+1 2 

9.37E-02  -3.50E-02  -9.37E-01  3.50E-01  2.78E+00 

1 8.20E+12 

-9.A8E-02  -2.50E-02  9.48E-01  2.50E-01  2.53E-01 

1 8.63E+12 

-1.64E-02  -2.50E-02  -5.49E-01  -8.36E-01  2.I5E+00 

1 2.02E+13 

1.64E-02  -2.50E-02  5.49E-01  -8.36E-01  9.90E-01 

I 9.90E+I2 

9.A8E-02  -2.50E-02  -9.68E-01  2.50E-0)  2.89E+00 

1 3.95E+12 

-9.89E-02  -I.50E-02  9.89E-0)  1.50E-01  I.5IE-01 

1 8.49E+12 

-2.59E-02  -1.50E-02  -8.65E-01  -5.02E-01  2.62E+00 

1 9.53E+12 

2.59E-02  -1.50E-02  8.65E-01  -5.02E-01  5.26E-01 

I 3.09E+12 

9.89E-02  -I.50E-02  -9.89E-01  1.50E-OI  2.99E+00 

I 2.S9E+1 2 

-9 . 99c -02  -5.00E-03  9.99E-01  5.00E-02  5.00E-02 

1 B.42E+12 

-2.95E-02  -5.00E-03  -9.S6E-01  -1.67E-01  2.97E+00 

I 8.32EH2 

2.95E-02  -5.00E-03  9.8AE-01  -1.47E-01  I.A8E-01 

1 2.30E+12 

9.99E-02  -5.00E-03  -9.99E-01  5.00E-02  3.09E+00 

1 2.20E412 


Figure  3.4  (continued)  MCCABE  output  for  the  coax  defined  in  Figure  3.2 


u.  , 

NMIX  = 4, 

PHI  = 0., 

XKT  = 0., 

JPRT  = 0, 

♦END 

UMIXTURE 

2A  = 2. 9E*01 , 0.,  0.,  0.,  0.,  0.,  0.,  0.f  0.,  0., 

UT  = 1.0,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

AT  = 6.354E+01,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

DENS  I T Y = 8.9, 

IZDM  = 29, 

♦ END 

1 ♦MIXTURE 

ZA  = 1.0,  6.0,  7.0,  8.0,  0.,  0.,  0.,  0.,  0.,  0., 

UT  = 9.41E-03,  5.6091 E-Ot , 1.3082E-01,  2.9886E-01,  0.,  0.,  0.,  0., 

+ 0. , 0. , 

AT  = 1.0,  1.2E+01,  1.4E+01,  1.6E+01,  0.,  0.,  0.,  0.,  0.,  0., 

DENSITY  = 1.42, 

IZDM  = 905, 

♦ END 

UMIXTURE 

ZA  = 5.0E+01 , 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

UT  = 1.0,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

AT  = 1 . 187E+02,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

DENSITY  = 7.2, 

IZDM  = 50, 

♦ END 

UMIXTURE 

ZA  = 4.7E+01 , 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

UT  = 1.0,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

AT  = 1.079E+02,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

DENSITY  = 1.05E+01, 

IZDM  = 47, 

♦ END 

Figure  3.5  (a)  Input  file,  TAPE5,  for  MCCABE  created  by  MCHELP  for  a 3-wire 
multi-conductor  cable. 
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1 


/ 


1 SATTPLN 
NPLANE  = 0, 

T - 0. , 0.,  0 > r 0. , 0 . , 0 . f 0.,  0.,  0.,  0.,  0.,  0.,  0 * , 0. , 0,,  0 

* , 0*1  0 * , 0. , 0 . , 

IPLANLB  = 0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 
SEND 

1SC0NDUCT 

RADC  = 1.8E-01,  5.0E-02,  5.0E-02,  5.0E-02,  0.,  0.,  0.,  0.,  0.,  0.,  0 
+ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

0., 

RADII  = 1 .50001 E-01 , 6.95E-02,  6.95E-02,  6.95E-02,  0.,  0.,  0.,  0.,  0. 
+ , 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

0. , 0. , 0. , 

RO  = 0.,  8.03E-02,  8.03E-02,  8.03E-02,  0.,  0.,  0. , 0.,  0.,  0.,  0., 
+ 0. , 0. , 0. , 0. , 0. , 0. , 0. , 0. , 0. , 

0., 

THET  = 0.,  6.0E+01,  1.8E+02,  3.0E+02,  0.,  0.,  0.,  0.,  0.,  0.f  0.,  0. 

+ , o.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

GAP  = 5.0E-03,  1.0E-08,  1.0E-08,  1.0E-08,  0.,  0.,  0.,  0.f  0,,  0.,  0 

a , 0*  j 0*|  Oaf  Oaf  Oaf  Oaf  Oaf  Oaf  Oaf 
Oa, 

NCON  = 4, 

TSHIELD  = 3.3E-02, 

HATLBLC  = 1,  1,  1,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0 
* * 

HATLBLD  = 2,  2,  2,  2,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0 

NY  =10, 

NATBAK  = C, 

NATLBLF  =3, 4, 4, 4, 

DFUSH=4*2.E-4, 

IPRT  = 0, 

SEND 


Figure  3.5  (b)  Input  file,  TAPE3,  for  MCCABE  created  by  MCHELP  for  a 3-wire 
multi-conductor  cable. 


1JNAHE1  = 5.0E-02,  5.0E-02,  5.0E-02,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0. 

t , Oa,  Oa,  Oa,  Oa,  Oaf  Oa,  Oa,  Oa, 

RH  = 8.03E-02,  8.03E-02,  8.03E-02,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

+ Oa,  Oa,  Oa,  Oa,  Oa,  Oa,  Oa,  Oa,  Oa, 

NHAX  = 3, 

UMAX  = 4, 

PER  = 3.5, 

THSLN  = 6!oE+01,  1.8E+02,  3.0E+02,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0. 

+ , Oa,  Oa,  Oa,  Oa,  O.f  Oa,  Oa,  Oa, 

IPRT  = 0, 

SEND 


Figure  3.5  (c)  Input  file,  TAPE15 , for  MCCABE  created  by  MCHELP  for  a 3-wire 
multiconductor  cable. 
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1  CIRCLES  8/3/7? 


MINIMUM  DIAGONAL  1.24996E-01 

CAPACITANCE  MATRIXIF/CH ) 

1 I 1.33434E-12 

I 2 5.80652E-13 

1 3 5. 80652E-1 3 

2 1 5. 80652E-1 3 

2 2 1.33434E-12 

2 3 5 .80652E- 1 3 

3 1 5 . 80652E- 1 3 

3 2 5 .80652E- 1 3 

3 3 1 . 33434E- 1 2 


1 DRIVPAK  8/3/77 


GEOMETRY  VARIABLES  TRANSFERRED  FROM  ATTNPAK 
NCON=  4 

RADC=  1.80E-01  5.00E-02  5.00E-02  5.00E-02 

GAP=  5.00E-03  1.00E-08  1.00E-08  1.00E-08 

RADD=  1 . 50E-01  6.95E-02  4.95E-02  6.95E-02 

RO=  0.  8.03E-02  8.03E-02  8.03E-02 

THET=  0.  6 . OOE+OI  1.80E+02  3.00E+02 

HATLBLC=  1111 

MATLBLD=  2222 

MATLBLF  = 3 4 4 4 

DFLASH=  2.00E-04  2.00E-04  2.00E-04  2.00E-04 

PHOTON  ENERGIES 
NE  = 1 

EPHOT=  5.00E+01 

TOTAL  EMITTED  CURRENT  FROM  SHIELD  AND  EACH  CONDUCTOR 
ISURF=  6.34E-08  1.38E-08  2.20E-08  1.38E-08 

TOTAL  INDUCED  CURRENT  ON  EACH  CONDUCTOR  DUE  TO  DEPOSITED  CHARGE 
I AREA=  -1.73E-08  -2.99E-08  -1.73E-08 

NORTON  DRIVER  TO  EACH  CONDUCTORISHORT  CIRCUIT  CURRENT)  (A/CM) 
ISC=  -3.44E-09  -7.89E-09  -3.44E-09 


Figure  3.6  MCCABE  output  for  the  three-wire  problem  defined  in  Figure  3.5. 
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APPENDIX  A 

REVIEW  OF  POTENTIAL  THEORY 


Consider  a set  of  N+l  parallel  conductors,  the  first  of  which  is  grounded,  separated 

by  a various  homogeneous  dielectrics.  A volume  charge  density  p(r)  is  deposited  in  the 

dielectric  at  the  points  r'.  The  potential  <p , the  solution  of  Poisson's  equation,  and  the 

3d)  3 G 

Green  function  G satisfy  the  following:  <p , G,  K->-  ~ and  K-*-  - — are  continuous,  G=0  on  the 

r on  r on 

conductor  boundaries,  and 

v2*  - P v2r  = - r')  . . 

’ * ’ ' E K+  ’ 5 G ~ E K+  ’ (A,1) 


where  e K-*  is  the  dielectric  constant  at  the  point  r.  From  Green's  theorem, 
o r 


/d2r  (' 


<})V2k-k;  - k-»gv2$ 

r r 


) - f £) 


we  obtain  a representation  for  4> , 


3G(r.r) 


/•  2 INTi  ~ a KJ  \ L , 

d r"p  (r')G(r'r)-e  V V.  /ds.  R*  — r— 1 

1 J 3 r.  dn. 


, (A. 3) 


where  V.  is  the  potential  of  the  i conductor.  The  normal  n.  at  the  surface  whose  in- 
J 3 

cremental  line  element  is  ds^ , is  directed  out  of  the  dielectric  volume.  A second  ap- 
plication of  Green’s  theorem  (with  G replacing  4>  in  Eq.  (A. 2))  establishes  a symmetry  of 

G, 


G(rr')  = G(r"r). 


The  induced  charge  on  conductor  i is 


Q = e fds . K-+ 
i o / i r.  3n. 

J it 

= Co/dsiKri  /d2r'p(?'>  jh 


I 77'  _ ...  - 


3 

3n . 
J 


K-+  G(r.r.) 
r . j 1 
J 


(A.  5) 


We  proceed  to  simplify  the  interpretation  of  Eq.  (A. 5).  First  we  define  the  capacitance 
matrix 


k.  . 


2 fds  3 i;-*.  fd 

° J i r LJ  j 


K>  G(r  r.). 

r.  3n . i i 

J J J 


(A. 6) 


Next,  considered  Eq.  (A. 3)  with  p = 0,  and  with  all  conductors  grounded  except  con- 
ductor i which  has  unit  potential.  The  potential  at  the  point  r',  which  we  now  call 
ip^(r"),  satisfies 


K+  i — G(r.r  ), 
r.  3n.  l o 

l l 


(A. 7) 


This  can  be  inserted  into  Eq.(A.5)  and  we  have 


h-f 


d2r ' 


p(r‘ 


)^(r') 


N 

E 

j=i 


'ij 


(A. 8) 


The  first  term  represents  the  induced  charge,  the  second  the  capacitive  charge,  for  the 
elemental  length  of  cable. 


EXPLICIT  EXPRESSIONS  INVOLVED  IN  LAPLACE  SOLUTION 

->  -* 

The  vectors  r = r,  <j>  and  p = p^9^  ha\e  the  shield  center  as  origin;  r 
trary  point  at  which  the  potential  is  evaluated,  p is  the  location  of  the 
center  with  respect  to  the  origin.  Define 

2 2 

r pi  2 

S = 2~  + aQ  - 2 rpt  cos  ($-0^, 

a 

o 


T = r2  + p±2  - 2rp1  cos  (<(>-6.^, 


Then 


U = arctan 


rsini|>  - p^  sin0^ 
rcosi(>  - p^  cos©^  ’ 


In  S/T,  n=0 

m"  t 


cos  (nU)" 
s in  (nU ) 


oo 


cos  ([ra-n^]  0^  - tn<J>) 
-sin  ([ra-n']  4> ^ - m4>) 


, n ^ 0 


Replacing  r p . and  p -*•  p ,-.,in  S,  T and  U above,  we  have  for  A . „ 
1 1 i in»i  n 

where  n=0,  n'=0 


2 2 
a - P, 

V.  ...  = 6...  In  -2 — + (1  - 64J  .)  . S 

in, i n ii  a^  ii  ' ln  — ; 


is  an  arbi- 
ith  conductor's 


where  n>0,  n'=0 
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U.S.  Army  Communications  Sys.  Agency 
ATTN:  CCM-AD-LB,  Library 

U.S.  Army  Foreign  Science  & Tech.  Ctr. 
ATTN:  DRXST- IS- 1 

DEPARTMENT  OF  THE  NAVY 

Naval  Research  Laboratory 

ATTN:  Code  7550,  J.  Davis 
ATTN:  Code  6701 

Naval  Surface  Weapons  Center 
ATTN:  Code  F31 

Office  of  Naval  Research 

ATTN:  Code  427,  H.  Mullaney 

Strategic  Systems  Project  Office 

Department  of  the  Navy 
ATTN:  NSP 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Geophysics  Laboratory 
ATTN:  C.  Pike 

Air  Force  Weapons  Laboratory 


ATTN: 

SUL 

ATTN: 

NT 

2 

cy 

ATTN: 

DYC 

2 

cy 

ATTN: 

NXS 

Deputy  Chief  of  Staff 
Research,  Development,  & Acq. 

Department  of  the  Air  Force 
ATTN:  AFRDQSM 

Rome  Air  Development  Center,  AFSC 
ATTN:  ESR,  E.  Burke 

Space  & Missile  Systems  Organization/DY 
Air  Force  Systems  Coimand 
ATTN:  DYS 

Space  A Missile  Systems  Organization/MN 
Air  Force  Systems  Command 
ATTN:  MNNG 
ATTN:  MNNH 
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DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Space  A Missile  Systems  Organization/SK 
Air  Force  Systems  Command 
ATTN:  SKF 

Strategic  Air  Command/ XPFS 
Department  of  the  Air  Force 

ATTN:  NRI-STINFO,  Library 
ATTN:  XPFS 


Hughes  Aircraft  Co. 

El  Segundo  Site 

ATTN:  MS  A1080,  W.  Scott 
ATTN:  MS  A620,  E.  Smith 

Institute  for  Defense  Analyses 
ATTN;  Classified  Library 


IRT  Corp. 


OTHER  GOVERNMENT  AGENCIES 

ATTN: 

Library 

ATTN: 

D.  Swift 

NASA 

Lewis  Research  Center 

JAVCOR 

ATTN 

N.  Stevens 

ATTN: 

E.  Wenaas 

ATTN 

C.  Purvis 

ATTN: 

Library 

ATTN 

Library 

JAYCOR 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

ATTN: 

R.  Sullivan 

Aerospace  Corp. 

Johns  Hopkins  University 

ATTN 

F.  Hai 

Applied  Physics  Lab. 

ATTN 

J.  Reinheimer 

ATTN: 

P.  Partridge 

ATTN 

V.  Oosephson 

ATTN 

Library 

Kaman  Sciences  Corp. 

ATTN: 

J.  Lubell 

Avco  Research  A Systems  Group 

ATTN: 

W.  Rich 

ATTN 

A830,  Library 

ATTN: 

Library 

Boeing  Co. 

ATTN:  P.  Geren 

University  of  California  at  San  Diego 
ATTN:  S.  De  Forest 

Computer  Sciences  Corp. 

ATTN:  A.  Schiff 


Dikewood  Industries,  Inc. 

ATTN:  Tech.  Lib. 

Dikewood  Industries,  Inc. 

ATTN:  K.  Lee 

EGAG  Washington  Analytical  Services  Center,  Inc. 
ATTN:  Library 

Eugene  P.  Deplomb 

ATTN:  E.  Deplomb 

Ford  Aerospace  A Communications  Corp. 

ATTN:  Technical  Library 
ATTN:  MS  G30,  D.  McMorrow 

General  Electric  Co. 

Space  Division 

ATTN:  VFSC,  Rm.  4230M,  J.  Peden 

General  Electric  Company- TEMPO 
Center  for  Advanced  Studies 
ATTN:  DAS  I AC 
ATTN:  W.  McNamara 

Hughes  Aircraft  Co. 

ATTN:  Technical  Library 


Lawrence  Livermore  Laboratory 
University  of  California 

ATTN:  Doc.  Con.  for  Technical  Information  Dept. 
Libra  ry 

Lockheed  Missiles  A Space  Co.,  Inc. 

ATTN:  Dept.  85-85 

Los  Alamos  Scientific  Laboratory 

ATTN:  Doc.  Con.  for  Reports  Library 

McDonnell  Douglas  Corp. 

ATTN:  S.  Schneider 

Mission  Research  Corp. 

ATTN:  R.  Stettner 
ATTN:  C.  Longmire 

Mission  Research  Corp. -San  Dieao 
ATTN:  Library 
ATTN:  V.  Van  Lint 

RAD  Associates 

ATTN:  L.  Schlessinger 

ATTN:  C.  MacDonald 

ATTN:  Technical  Information  Center 

Rockwell  International  Corp. 

ATTN:  Library 

Sandia  Laboratories 
Livermore  Laboratory 

ATTN:  Doc.  Con.  for  T.  Dell  in 

Sandia  Laboratories 

ATTN:  Doc.  Con.  for  3141 


